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Abstract 

We have developed a near-infrared camera called ANIR (Atacama Near-InfraRed camera) for the University of 
Tokyo Atacama Observatory 1.0-m telescope (miniTAO) installed at the summit of Cerro Chajnantor (5,640 m above 
sea level) in northern Chile. The camera provides a field of view of 5(1 x 5(1 with a spatial resolution of 0"298 
pixel 1 in the wavelength range from 0.95 to 2.4 //m, using Offner relay optics and a PACE HAWAII-2 focal plane 
array. Taking advantage of the dry site, the camera is capable of hydrogen Paschcn-o (Pact, A = 1.8751 pm in air) 
narrow-band imaging observations, at which wavelength ground-based observations have been quite difficult due to 
deep atmospheric absorption mainly from water vapor. We have been successfully obtaining Pact images of Galactic 
objects and nearby galaxies since the first-light observation in 2009 with ANIR. The throughputs at the narrow-band 
filters (A' 1875, A191) including the atmospheric absorption show larger dispersion (~ 10%) than those at broad¬ 
band filters (a few percent), indicating that they are affected by temporal fluctuations in Precipitable Water Vapor 
(PWV) above the site. We evaluate the PWV content via the atmospheric transmittance at the narrow-band filters, 
and derive the median and the dispersion of the distribution of the PWV of 0.40 ± 0.30 and 0.37 ± 0.21 mm for 
the A' 1875 and A191 data, respectively, which are remarkably smaller (49 ± 38% for A' 1875 and 59 ± 26% for 
A191) than radiometry measurements at the base of Cerro Chajnantor (an altitude of 5,100 m). The decrease in 
PWV can be explained by the altitude of the site when we assume that the vertical distribution of the water vapor 
is approximated at an exponential profile with scale height within 0.3-1.9 km (previously observed values at night). 

We thus conclude that miniTAO/ANIR at the summit of Cerro Chajnantor indeed provides us an excellent capability 
for a ground-based Pao observation. 

Key words: site testing — infrared: Paschcn-a — infrared: general — instrumentation: detectors 


1. Introduction 

1.1. Paa emission as a tool for Unveiling the Dust-obscured 
Universe 

To explore the formation and the evolution of the present- 
day “normal” galaxies like our Galaxy, it is important to char¬ 


acterize where and how star formation occurs in a galaxy. In 
the local universe, it is well known that star formation rate 
(SFR) correlates with dust opacity (Takeuchi et al. 2010; 
Bothwell et al. 2011), where intensely star-forming regions 
and galaxies are obscured by a huge amount of dust, and be¬ 
come optically thick (Alonso-Herrero et al. 2006b; Piqueras 
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Fig. 1. Atmospheric transmittances in the NIR wavelength range at 
TAO {top), Mauna Kea {middle, 4,200 m above sea level) and Paranal 
{bottom, 2,600 m), simulated using an atmospheric transmission model, 
ATRAN (Lord 1992) with a spectral resolution of ~ 50000. Median 
PWV of 0.5 (Giovanelli et al. 2001b), 1.8 (Otarola et al. 2010), and 
2.4 mm (Kerber et al. 2010) are assumed for TAO, Mauna Kea, and 
Paranal, respectively. Remarkable improvements of the transmittance 
are seen at around 1.9 pm as well as at gaps between conventional 
broad atmospheric windows such as J (1.28 pm), H (1.67 pm), and 
TCs-band (2.15 pm). 


Lopez et al. 2013). Therefore, the UV continuum emission 
and hydrogen recombination lines in the optical ranges such as 
Ha at 0.6563 p m and H/l at 0.4861 pm emitted from those 
dusty star-forming regions are easily attenuated, and then it 
would be complicated to correctly understand the distribution 
of star-forming regions in the whole galaxy, especially when 
the galaxy has a patchy distribution of dust (Garcfa-Marm et al. 
2006). 

On the other hand, Paschen-a (Paa) at 1.8571 pm (in air; 
1.8756 /rm in vacuum), the strongest emission line in the near- 
infrared (NIR, A ~ 0.8-2.5 p m) wavelength range, is less af¬ 
fected by dust thanks to its longer wavelength than the opti¬ 
cal lines (Kennicutt 1998), and becomes a powerful and di¬ 
rect indicator of SFR, especially in dusty regions. In particu¬ 
lar, the observed Paa emission becomes stronger than Ha at 
E{B - V) > 1.2 (A v > 3.7) and than Br 7 at E(B - V) < 
28.0 (Ay < 86.2) with an assumption of a Milky-Way-like ex¬ 
tinction curve and A V /E(B — V) = 3.08 (Pei 1992), while 
the intrinsic Paa luminosity is 0.12 of Ha, 2.14 of Pa/3 (1.2818 
/rm in air), and 12.4 of Br 7 (2.1655 /rm in air) for Case B re¬ 
combination with an electron temperature of 10 4 K and density 
of 10 4 cm -3 (Osterbrock 1989). 

However, Paa observations from the ground have been quite 
difficult so far due to deep atmospheric absorption features 
around their wavelength, mostly caused by water vapor, and 
hence most Paa observational studies have been limited to 
those by Near Infrared Camera and Multi-Object Spectrometer 
(NICMOS, Thompson et al. 1998) on board the Hubble Space 


Table 1 . Specification of the ANIR NIR Unit. 


Wavelength 

0.95-2.4 pm 

Detector 

PACE HAWAII-2 

Pixel format 

1024 x 1024* 

Pixel pitch 

18.5 pm 

Field of view 

5 A x 5.'1 

Pixel scale 

0"298 pixel’ 1 

Broad-band filters ^ 

Y, J, H, K s 

Narrow-band filters 

AT28, (V1875, iV 191, 7V207 


* The actual format size of HAWAII-2 FPA is 2048 X 2048 pixels, of 
which we use only 1024 X 1024 pixels. 

' The J, H, and K a filters are those of Mauna Kea Observatories 
(MKO) infrared filter set (Simons & Tokunaga 2002; Tokunaga et al. 
2002 ). 

Telescope ( HST ). For example, Alonso-Herrero et al. (2006b) 
use the NICMOS Paa narrow-band (F187N and F190N) imag¬ 
ing data to study star formation properties of dusty star-forming 
galaxies in the local universe with a high spatial resolution, and 
find the compact distributions of Paa emission along various 
galaxy structures such as nucleus and spiral arm. They also es¬ 
tablish a linear empirical relation between Paa and 24 /im lu¬ 
minosity as well as the total infrared luminosity as an indicator 
of SFR for those galaxies. In addition, a wide-held Paa imag¬ 
ing survey of the Galactic Center (Wang et al. 2010; Dong 
et al. 2011) produces a high spatial-resolution map of stars and 
an ionized diffuse gas with a possible new class of massive stars 
not associated with any known star clusters. Recently, a NIR 
camera and spectrograph, FLITECAM (Mclean et al. 2012), 
for NASA’s Stratospheric Observatory for Infrared Astronomy 
(SOFIA) has been developed, and will support Paa imaging 
observations with two narrow-band filters centered at 1.875 and 
1.900 /an. 

As described above, Paa observations are quite useful for 
a variety of targets obscured by dust and can provide a new 
insight into those studies, yet the amount of observing time 
and facilities available for that purpose are still very limited. 
Therefore, we have developed a new instrument at an appro¬ 
priate astronomical site for ground-based Paa observations. 

1.2. ANIR: Ground-based Paa Imager for the miniTAO 1.0-m 
Telescope 

The University of Tokyo Atacama Observatory (TAO, Yoshii 
et al. 2010) 1.0-m telescope, called miniTAO, is an opti¬ 
cal/infrared telescope installed in 2009, located at the summit 
of Cerro Chajnantor (an altitude of 5,640 m or 18,500 ft above 
sea level) in northern Chile as the world’s highest astronomical 
observatory (Sako et al. 2008a; Minezaki et al. 2010). The 
site was expected to be very dry from satellite data reported by 
Erasmus & Sarazin (2002) (precipitable water vapor, PWV, of 
0.5 mm at 25th percentile) and from radiosonde measurements 
by Giovanelli et al. (2001b) (0.5 mm at a median above an 
altitude of 5,750 m). Figure 1 shows a comparison of simu¬ 
lated atmospheric transmittances between the TAO site and the 
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other sites at lower altitudes, in which there are remarkable im¬ 
provements of the transmittance at gaps between conventional 
broad atmospheric windows, especially around 1.9 pm, i.e., 
Pao wavelength. 

The Atacama NIR camera (ANIR) 1 is one of the instruments 
for the Cassegrain focus of the miniTAO telescope. It is capa¬ 
ble of wide-held Paa narrow-band (NB) imaging observations 
as well as ordinary broad-band imaging one. It also has a ca¬ 
pability of optical imaging and slitless spectroscopic observa¬ 
tions simultaneously with the NIR imaging using a retractable 
dichroic mirror. We have achieved first-light observations us¬ 
ing NB filters targeted on the Paa emission line (Figure 4) in 
2009, and now ANIR is in operation mainly for Paa observa¬ 
tions of the Galactic center/plane and nearby starburst galaxies 
in order to understand the star formation activities hidden in 
thick dust clouds (e.g., Tateuchi et al. 2012b; Komugi et al. 
2012 ). 


In this paper, we describe the overall design of ANIR in 
Section 2, its performance evaluated through laboratory tests 
and actual observations at the TAO site in Section 3, and a site 
evaluation through the Paa imaging data in Section 4. We re¬ 
fer the reader to Tateuchi et al. (2012a) and Tateuchi et al. 
(2014) for a detailed description of data reduction processes 
and quantitative analysis of Paa data. 

2. Instrument 

2.1. Near-infrared Imager Unit 

Here we describe the design of the NIR unit providing broad¬ 
band and NB imaging functions in the NIR. A hardware design 
of the optical unit is described in Section 2.2. 

2.1.1. Overview 

Figure 2 shows a cross-sectional view of ANIR, and a block 
diagram of the hardware is shown in Figure 3. The specifi¬ 
cations of the NIR unit are summarized in Table 1. The unit 
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ANIR web site: http://www.ioa.s.u-tokyo.ac.jp/kibans/anir_en/ 























































































































































































4 


M. Konishi et al. 


[Vol., 



Fig. 3. Hardware block diagram of ANIR. 


covers a field of view (FoV) of 5'1 x 5'1 with a spatial reso¬ 
lution of 0"298 pixel -1 using an engineering grade Producible 
Alternative to CdTe for Epitaxy (PACE) HAWAII-2 array de¬ 
tector, which is a 2048 x 2048 pixel HgCdTe NIR focal plane 
array (FPA) manufactured by Teledyne Scientific & Imaging 
LLC. Note that only a single quadrant with 1024 x 1024 pix¬ 
els is used. The filter set consists of four standard broad-band 
filters (Y, J, H, and K s ) and four NB filters (A 128, A1875, 
TV 191, and A207) as listed in Table 2. The A1875 filter has 
a central wavelength of A c = 1.8759 pm and a bandwidth of 
AA = 0.0079 pm, which covers the Paa (A =1.8751 pm) with 
radial velocities within ~ —580 to +680 km s _1 . On the other 
hand, the TV191 filter (A c = 1.911 pm, AA = 0.033 pm) corre¬ 
sponds to redshifted Pact lines with recession velocities of cz 
~ 2900-8200 km s _1 . The A191 filter is also used for taking 
off-band (continuum) data for the TV 1875 data, and vice versa. 

Figure 4 shows details of the atmospheric transmittances 
shown in the top panel of Figure 1 at the wavelength range 
of interest for the Paa NB filters. While there are absorption 
features mostly by water vapor even at the altitude of 5,640 m, 
several windows with a high transmittance exist. The average 
atmospheric transmittances at the PWV of 0.5 mm are approx¬ 
imately 0.48 and 0.64 within the bandpasses of the A1875 and 
A191 filters, respectively. 


2.1.2. Cryogenics 

The cryostat is a compact cube with approximately 260 mm 
on a side, in which Offner relay optics (Section 2.1.3), a fil¬ 
ter box with two wheels, and the HAWAII-2 FPA are housed, 
as shown in Figure 2. All the components in the cryostat are 
cooled down to 70 K to reduce thermal radiation, especially 
for observations in the A s -band. A single-stage closed-cycle 
mechanical cooler with a cooling capacity of 25 W at 77 K is 
equipped. A cold head anti-vibration mount and bellows are 
inserted between the cold head and the cryostat to reduce vi¬ 
bration caused by the cold head. In addition, a Teflon plate and 
a sapphire plate are inserted to electrically insulate the cryostat 
from the cold head. It takes approximately 24 hours to cool 
down and stabilize the cryostat (AT < 0.05 K at the detector 
box) from the ambient temperature down to 70 K. 

The two filter wheels are installed just after the focal plane of 
the telescope. Each filter wheel has five slots with tp 34 mm. A 
stepping motor controls its rotation, and neodymium magnets 
and a hall-effect sensor sense the position of the slots. 

2.1.3. Optics 

Reflective Offner relay optics are employed for re-imaging, 
consisting of two (concave primary and convex secondary) 
spherical mirrors. The primary mirror forms an image of the 
telescope pupil on the secondary mirror, which works as a cold 
Lyot stop. The specifications of the mirrors are summarized in 
Table 3. These mirrors are gold-coated to achieve high reflec- 
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Table 2. NIR narrow-band filters. 


Name Wavelength [/itn] Targeted Line 

Specification Measurement* 



Center (A c 1 ') 

Width (AA*) 

Center 

Width 


AT128 

1.2818 

0.0210 

1.2814 

0.0217 

Pa/? (1.2818 pm) 

ZV1875 

1.8751 

0.0080 

1.8759 

0.0079 

Pact (1.8751 pm) 

W191 

1.9010 

0.0306 

1.9105 

0.0329 

Pact oH -hancI ! 

/V207 

2.0750 

0.0400 

2.0742 

0.0391 

Civ (2.078 pm) 


* Measured at 77 K. 

T Isophotal wavelength (Tokunaga & Vacca 2005). 

I Full-width at half maximum. 

§ Also corresponding to Pact rest-frame wavelength emitted from extra-galactic objects within c z ~ 
2900-8200 km s -1 . 



1.86 1.87 1.88 1.89 

Wavelength [//in] 



1.88 1.90 1.92 1.94 

Wavelength [pm] 


Fig. 4. Atmospheric transmittances at the TAO site simulated with the PWV of 0.5 mm are shown for the 771875 (left) and N 191 (right) filters, re¬ 
spectively. The shaded region in the top panels represents the effective transmittance including the atmospheric transmittance (solid line) and the filter 
transmittance measured at 77 K (shown in the bottom panel). The vertical dashed line indicates the position of Pact wavelength (1.8751 pm). The average 
atmospheric transmittance within the bandpass (AA) around the center (A c ) listed in Table 2 is 0.48 and 0.64 at the .V1875 and .Y191 filters, respectively. 
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Fig. 5. Spot diagrams at A = 1.65 /rm without (left) and with (right) the dichroic mirror. Diagrams at the center of the FoV and positions at a 5 X 5 grid 
with a step of ±1.'32 from the center are shown. A side of each box corresponds to 90 / /m or 1”5. The circle corresponds to the Airy disk at A = 1.65 
/im. It is worth mentioning that no significant image degradation caused by the insertion of the dichroic mirror is seen throughout wavelengths between 
0.95 and 2.4 /im. 


Table 3. Specification of the Offner relay optics. 


Primary mirror 

Radius of curvature 

140 mm 


Effective diameter 

90 mm 

Secondary mirror 

Radius of curvature 

70 mm 


Effective diameter 

9 mm 

Offset of optical axis 


24 mm 


tivity. 

Figure 5 shows spot diagrams at A = 1.65 ft m across the 
FoV. Sharp image quality is achieved at any positions within 
the FoV. As shown in the right-hand panel of Figure 5, even 
with a dichroic mirror inserted for simultaneous optical-NIR 
observations (Section 2.2), spot sizes are still smaller than or 
comparable to the size of the Airy disk. Figure 6 shows encir¬ 
cled energy distributions, the fraction of energy in a given ra¬ 
dius to the total energy from a point source, for the center and 
off-center spots with the dichroic mirror inserted, compared to 
that of a diffraction-limited spot. We confirm that the diameter 
encircling the 80% energy is less than 1.7 pixels (~ 31.5 /fin) 
or 0"5 across the FoV. 

2.1.4. Data Acquisition and Control System 

ANIR is operated by two Linux PCs (named uni and 
uni2)\ one is dedicated to control HAWAII-2 FPA (TAO Array 
Controller, TAC), and the other handles all the other tasks, in¬ 
cluding operation of the filter wheels and the optical unit (see 
Section 2.2), acquisition of house keeping information such as 
temperature and vacuum pressure, and management of various 



Radius from centroid |//m 


Fig. 6. Encircled energy distributions of the point spread functions of 
the NIR unit at the positions (O.'O, —2'6) and (O.'O, +2'6) from the 
center of the FoV (O.'O, O'O) (thick lines) in the case that the dichroic 
mirror is inserted, compared to that of a diffraction-limited image (thin 
solid line). 
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Obs. Commands 


ANIR 


Observatory Server 


MySQL Database 


Telescope/Dome Status 
Weather Status 
Object Information 


; FITS 
Tel. Commands 1 Generation 


Latest [}B data 


Tel. Commands 


Telescope Server 


Linux (uni) 


> ^ \ 


\ 


Linux (uni2) 

Inst. Commands NIR Filter Operation 

_ ^ Optical Observation 

^_ Instrument Status Monitor 

Web Server 


Latest DB data 


Clock 


1024 x 1024 
pixel Data 


CCD Commands 

-> f 


MySQL Database 


Instrument 
Temperatures/Vacuum 
Opt/NIR Filter Status 
HAWAII-2 Parameters 
Dichroic Mirror Status 


1056 x 1027 
pixel Data 


Fig. 7. Software block diagram of ANIR. 


instrument status. The TAC system (see Sako et al. 2008b for 
more details) is a high-speed and flexible array controller using 
a real-time operating system instead of conventional dedicated 
processor devices such as digital signal processors (DSPs). 
One core of a dual-core CPU in uni is assigned to the real¬ 
time data processing, and engages in clock pattern generation 
and processing of acquired frame data. This allows us to con¬ 
trol the FPA in real time without affecting the performance of 
the operating system and other software. A rate to read a pixel 
(Pixel rate) can be set to 3 to 8 /is. Considering increasing 
readout noise with faster readout (slower Pixel rate), the rate of 
4 iis is usually used, corresponding to a readout time of 1024 
x 1024 pixels of 4.2 s, which is the shortest exposure time of 
the NIR unit. 

Figure 7 summarizes the functions and gives a com¬ 
mand/data flow diagram of the control system. To efficiently 
handle status data increasing with time, we adopt a relational 
database management system, MySQL. An successful exam¬ 
ple of the application of MySQL in astronomical instrumenta¬ 
tion is described in Yoshikawa et al. (2006). Several database 
tables are prepared in MySQL, and information concerning in¬ 
strument status (filters used, insertion of the dichroic mirror, 
temperatures/pressures) is stored into the respective tables ev¬ 
ery minute. When the exposure command is executed, a FITS 
(Flexible Image Transport System) header is constructed by 
collecting the latest information from those tables. 

2.2. Optical Unit 

To benefit from the advantages of the site, such as good see¬ 
ing condition with 0"7 at L-band at a median (Motohara et al. 
2008; see also Giovanelli et al. 2001a), ANIR is capable of 
optical and NIR simultaneous imaging observations by insert¬ 
ing a dichroic mirror in front of the entrance window of the 


cryostat of the NIR unit (the upper part of Figure 2). 

The specifications of the optical unit are summarized in 
Table 4. The unit covers a FoV of 6.'0 x 5(9 with a spa¬ 
tial resolution of 0."343 pixel -1 using a commercial, peltier- 
cooling CCD camera unit incorporating a E2V backside- 
illuminated CCD, Proline PL4710-1-MB, fabricated by Finger 
Lakes Instrumentation. Four Johnson-system broad-band fil¬ 
ters ( B , V, R, and I) for imaging and a low-resolution trans¬ 
mission grating (grism) with 75 lines mm -1 and a blaze angle 
of 4.3° for slitless spectroscopy are available. The dichroic 
mirror has a dimension of 50 x 71 x 10 mm with a wedge 
angle of 0.78° to minimize astigmatism in the NIR image. Its 
reflectance at 0.4-0.86 /itn is higher than 0.9 with an incident 
angle of 45°, while its transmittance at 0.95-2.4 pm is kept 
higher than 0.9. Figure 8 shows spot diagrams of the optical 
unit at the final focal plane. Due to chromatic aberration of the 
re-imaging refractive optics, the focal position is offset by ~ 
0.2 mm in the /-band, which can be effectively compensated 
as shown in the right-hand panel of Figure 8 by shifting the 
camera optics mounted on a linear stage. The dichroic mirror is 
mounted on a linear stage, and can be retracted to derive better- 
quality (higher efficiency, lower thermal background) NIR data 
when no simultaneous optical observation is required. 

3. Performance 

In this section, we describe the imaging performance of the 
NIR and optical units obtained through observations carried out 
in 2009-2011. 

3.1. Detector Performances 

We measure the conversion factor of the HAWAII-2 FPA 
readout system using flat frames with increasing exposure time 
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Fig. 8. Left : Spot diagrams of the optical unit at the focus position. For each panel, a diagram for six wavelengths (0.44, 0.53, 0.62, 0.71, 0.80, and 0.90 
/xm) and four field positions (the center, and 2'04, 2'88, 4.'08 from the center) is shown. A side of each box corresponds to 100 /im or 2. // 5. Right: Same 
as the left-hand panel, but with the focus position offset by 0.2 mm to improve the performance in the I -band (0.81 /xm). 


in the laboratory, and obtain 3.2 tC ADU -1 from the calcu¬ 
lated mean and dispersion of pixel values. The linearity of the 
FPA is evaluated from other flat images taken by increasing ex¬ 
posure time. We ensure that the linearity is kept within ~ 1% 
up to 18,000 ADU or 57,600 e . The readout noise is mea¬ 
sured on the telescope in dark frames with the shortest exposure 
(i.e., 4.2 s) which are images with the J-band and V1875 filter 
inserted, where no photons fall on the FPA because those band- 
passes do not overlap each other. We obtain 16.6 e - as the rms 
of their median values with the Correlated Double Sampling 
(CDS) readout method. A multiple readout (“up-the-ramp" 
sampling) method is also examined, and we obtain 13.9 e - 
with 16 readouts. The dark current derived from dark frames 
with 300-500 s exposures is 0.28 e - s -1 pixel -1 . We should 
mention that the dark current of F1AWAII and F1AWAII-2 FPAs 
have persistence and that the level of the persistence tends to be 
dependent on incident flux (Hodapp et al. 1996; Finger et al. 
1998; Motohara et al. 2002). For example, the persistence 
of a dark current pattern would appear strongly on frames with 
a NB filter taken after exposures of high background, such as 
ones with a broad-band filter. The TAC system continues to 
reset the FPA with an interval of ~ 0.1 s in order to reduce the 
persistence effect even if no exposure command is issued. 


Table 4. Specification of the ANIR Optical Unit. 


Wavelength 

4000-8600 A 

CCD Unit 

FLI Proline PL4710-1-MB 


(E2V backside-illuminated CCD) 

Pixel format 

1056 x 1027 

Pixel pitch 

13.5 pm 

Field of view 

6.'0 x 5(9 

Pixel scale 

0"343 pixel -1 

Filters 

B, V, R, I (Johnson System) 

Grism 

75 lines mm -1 , blaze angle = 4.3° 


The performance of the optical CCD are also evaluated in 
a similar manner to that described above. We take flat im¬ 
ages with increasing exposure time in the laboratory, and obtain 
the conversion factor of 1.818 e~ ADU -1 from the calculated 
mean and dispersion of pixel values. The readout noise is mea¬ 
sured from bias frames to be 14.5 e - rms. We obtain the dark 
current of 0.011 e - s -1 pixel -1 at the CCD temperature of ~ 
220 K. 

3.2. Imaging Quality 

We evaluate the image quality of the NIR unit using bright 
(but unsaturated) and isolated field stars in images taken under 
good seeing condition. By matching positions on the image 
( Xi , Y,) of about 20 point sources located over the FoV with 
coordinates (ctj, 6,) obtained from the Two Micron All Sky 
Survey (2MASS) Point Source Catalog (PSC, Skrutskie et al. 
2006), we obtain a pixel scale of 0"298 pixel -1 with a fitting 
uncertainty of ~ 0." 1 pixel -1 . Image distortion is confirmed 
to be negligible (< 1 pixel) over the FoV, as designed. We 
derive a typical FWHM size of 2.41 ± 0.24 pixels or Of'72 ± 
0"07 in the K s -band where ellipticities of the stars are small 
and uniform (0.08 ± 0.04). 

The image quality of the optical unit is evaluated in a similar 
manner. Raw images have non-negligible image distortion, es¬ 
pecially at their corners. We evaluate and correct for the distor¬ 
tion in the similar manner for the NIR unit as described above, 
by matching the positions of stars with those in the HST Guide 
Star Catalog (GSC, Lasker et al. 2008) version 2.3 and fitting 
them to the latter using fifth-order polynomials to take into ac¬ 
count the image distortion. The mapping uncertainties are ~ 
0.3 pixels. We obtain a pixel scale of 0"343 pixel -1 at the 
center of the FoV for all the broad-band filters (B , V, R, and 
I), and a typical FWHM size of 3.04 ± 0.22 pixels or l."04 ± 
0"08 with a homogeneous ellipticity distribution (0.04 ± 0.08) 
in the F'-band. 

Figure 9 shows the seeing distribution for images taken in 
the latter half of the observation in 2010 and the former half in 
2011. The seeing size (FWHM) is measured using the IRAF 
psfmeasure task. The median seeing size of ~ 1" or less is 
obtained for all the bands. When we consider the diffraction- 
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Fig. 9. Seeing statistics in the optical ( left ) and NIR ( right ) for miniTAO/ANIR. The histograms are normalized at their peak, and their median value is 
represented as the vertical dashed line. A number below the respective filter name indicates the number of frames used for the statistics. 


limited size and Hartmann constant (0"19) of the telescope op¬ 
tics (Minezaki et al. 2010) and suppose that some data may be 
undersampled (FWHM < 2 pixel) due to good seeing condi¬ 
tion, the actual seeing distributions at the TAO site are expected 
to have a peak at a smaller size. 

3.3. Throughputs and Limiting Magnitudes 

Throughput is defined as the ratio of the number of photons 
coming from an object into the Earth atmosphere of an aper¬ 
ture of the telescope to that of electrons detected, where atmo¬ 
spheric extinction is included. The former is calculated with 
the effective collecting area of the telescope (~ 0.78 m 2 ) and 
the bandwidth of a filter used. The latter is simply obtained by 
multiplying a count rate of an object (ADU s -1 ) by the detector 
conversion factor. 

We use optical standard stars (Landolt 1992) for the optical 
unit and 2MASS PSC stars within 11-17 mag and with a me¬ 
dian photometric error of ~ 0.1 mag for the NIR unit for the 
measurements. The results are shown in Figure 10 and sum¬ 
marized in Table 5. Since the 2MASS catalog has only the 
J, H, and /\ s -band data, we calculate the magnitudes of the 
PSC stars for the Y -band and NB filters as follows. For the 
Y -hand, we linearly (in logarithmic scale) extrapolate the J- 
band (1.275 pm) magnitude to the F-band (1.038 pm) with a 
slope between the J and H -bands (1.673 pm) for each star. 
For the TV 128 filter, we use the J-band magnitude. The lower 
throughput of the TV 128 filter compared to the J-band is at¬ 
tributed to the combination use of the J-band filter as a ther¬ 
mal blocker for observations with the iV128 filter. We calcu¬ 
late magnitudes at the TV1875, iV191 and TV207 filters by lin¬ 
early (in logarithmic scale) interpolating their H and A' s -band 
(2.149 pm) magnitudes. We examine the influence of stellar 
photospheric features within the wavelength ranges of inter¬ 


est on the interpolation. At first, we use J — H and H — K s 
color-color diagram to explore which kind of stars dominate 
the 2MASS PSC, and derive median colors of J — H ~ 0.7 
and H — K s ~ 0.2, which correspond to intrinsic colors of K- 
to M-type stars (Bessell & Brett 1988; Stead & Hoare 2011) 
or would be those of G- to K-type stars with dust extinctions 
of Ay ~ 1-2 mag for data in low Galactic latitudes such as 
Galactic plane and Large/Small Magellanic Clouds (Imara & 
Blitz 2007; Dobashi et al. 2009). Then, we evaluate the 
influence of spectral features of those medium-mass stars by 
using a stellar atmosphere model of Kurucz (1993) containing 
O- to M-type stars. We find no significant differences (< 0.01 
mag) between magnitudes derived by interpolation and those 
derived by direct integration of the model spectrum convolved 
with the transmission curve of the NB filter, although those 
stars have hydrogen absorption lines, of course including Paa. 
Therefore, the magnitudes at those NB filters we use hereafter 
are thought to be reliable enough for the performance verifica¬ 
tion of those NB filters and the site evaluation (Section 4). The 
resulting throughputs of the TV1875 and iV 191 filters are sig¬ 
nificantly lower than those of the H and AT s -bands, and have 
larger dispersion by a factor of two, which are thought to be 
due to the lower atmospheric transmittance combined with the 
airmasses of the observations and temporal fluctuations in the 
atmospheric absorption mainly from water vapor (see Figure 4 
for the transmittance with a PWV of 0.5 mm at an airmass of 
1.0). We evaluate the PWV content above the TAO site by us¬ 
ing the throughputs of those filters in Section 4. The through¬ 
put of the TV207 filter is consistent with that of the /C-band 
when taking into account the difference in their average filter 
transmittances: ~ 78% (TV207) and ~ 92% ( K s ). It should 
be noted that, although there is several night-sky OH emission 
lines within the bandpasses of the NB filters (e.g., Rousselot 
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Table 5. Throughputs and limiting magnitudes*. 



Throughput 

[%] 

Limiting magnitude [mag]* 

@ 60 s @ 600 s 

Sky brightness 
[mag arcsec” 2 ] 

Tblip* 

[s] 

B 

17.5 

21.3 

23.5 

22.0 

2600 

V 

29.2 

21.4 

23.2 

20.8 

800 

R 

29.4 

21.5 

23.2 

20.3 

450 

I 

20.1 

20.6 

22.0 

18.1 

70 

Y 

14.3 

18.8 

20.1 

16.8 

130 

J 

18.9 

19.5 

20.9 

16.4 

50 

H 

29.1 

18.4 

19.7 

14.6 

10 

K s 

30.3 

18.7 

20.0 

15.3 

20 

iV128 

17.0 

16.9 

18.3 

14.7 

80 

TV 1875 

18.0 

15.6 

16.9 

13.5 

100 

iV191 

17.7 

16.9 

18.3 

15.0 

120 

N201 

27.7 

17.6 

18.9 

15.5 

100 


* All magnitudes are described in the AB magnitude system. 

' For a point source with S/N = 5 and aperture. 

1 Approximate exposure time required for background-limited performance (BLIP). We define BLIP as 
when the background noise becomes twice as large as the readout noise. 



Fig. 10. Throughputs of the ANIR optical and NIR unit on the mini- 
TAO telescope. The filled (open) circles show those of the broad-band 
(narrow-band) filters. The horizontal error bars represent the band- 
widths, and the vertical error bars show la dispersions. Note that at¬ 
mospheric extinction is not corrected. 

et al. 2000), we consider their influence on our result to be 
small because we perform aperture photometry with local sky 
subtraction for each source. 

The limiting magnitudes for the optical and NIR units are 


estimated using the dark current, readout noise, throughputs 
obtained above, and sky brightness (photons cm -2 s -1 /jm _1 
arcsec -2 ). The sky brightness is measured on a “sky” frame 
made by combining unregistered subsequent frames with ob¬ 
jects masked. Table 5 gives the estimated limiting magnitudes 
for a point source with S/N = 5, (j> 1"5 aperture, and two dif¬ 
ferent exposure times (60 and 600 s). Also listed is an approx¬ 
imate exposure time (Tblip) required to achieve background- 
limited performance (BLIP) that means the S /N of data is dom¬ 
inated by the background noise. We define BLIP as when the 
Poisson noise of the sky background becomes twice as large as 
the readout noise. 

3.4. Paa Imaging Performance 

Here we briefly demonstrate the performance of the mini- 
TAO/ANIR ground-based Paa observations. The reduced Paa 
emission image of the nearby starburst galaxy, IC5179 taken 
with the N191 filter, is shown in Figure 11, demonstrating that 
the TAO site indeed provides access to the Paa emission line 
from the ground. Star-forming regions associated with the spi¬ 
ral arm and the nucleus are clearly seen. For a comparison 
of its sensitivity, the same galaxy taken by the //5T/NICMOS 
F190N filter is also shown in the right-hand panel of Figure 
11, from which the F187N image is subtracted as a continuum 
emission. A quantitative comparison confirms the consistency 
of Paa fluxes obtained with ANIR and NICMOS within an ac¬ 
curacy of 10% (Tateuchi et al. 2012a). We refer the reader 
to Tateuchi et al. (2012a) and Tateuchi et al. (2014) for 
the data reduction processes, verification, quantitative analy¬ 
ses, and comparison of the performance with NICMOS, of Paa 
emission line data taken by the N 1875 and 191 filters. 
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Fig. 11. Postage stamps of the nearby starburst galaxy, IC5179 (a = 22 h 16 m 09? 10, <5 = —36 d 50 m 37F4, J2000.0), as an example demonstrating the 
ANIR Pan imaging performance. The left panel shows a 1(5 X 1.'5 composite image (R: continuum-subtracted Pact, G: /C-band, B: 77-band), oriented 
with north up and east to the left. The middle panel shows the continuum-subtracted Pact image (on-source exposure time of 4200 s). For comparison, 
a Pact emission line image taken by //.S77NICMOS is also shown in the right panel (exposure time of 960 s) from which the F187N image is subtracted 
from the F190N image as a continuum emission. Note that the NICMOS image covers only approximately the central 0.'4 X 0(4 region of the galaxy. 


4. Site Evaluation through Part Narrow-band 

Observations 

Finally, we evaluate how the site is suitable for the infrared 
astronomy in terms of PWV by using the N 1875 and A' 191 
data taken between 2009 and 2011. We note that our approach 
using the NB filters around 1.9 /im which are sensitive to PWV 
is independent of any previous studies carried out at (the vicin¬ 
ity of) the Chajnantor Plateau (e.g., Giovanelli et al. 2001b; 
Erasmus & Sarazin 2002; Peterson et al. 2003; Tamura et al. 
2011 ). 

The throughputs derived in Section 3.3 (hereafter ?/ 0 | JS ) con¬ 
sist of the following factors: (i) the reflectivity of the telescope 
mirrors (7V,,|), (ii) the system efficiency of ANIR excluding 
the filter transmittance (Tanir), (iii) the filter transmittance 
(^filter), (iv) atmospheric transmittance (T^m V ) at the zenith 
which is dependent on PWV. Then, the throughput is described 
as: 

band rr iband v /Tiband v Riband v //tiPWV, band\X xi\ 

Vobs — -*Tel x 1 ANIR x 1 filter x V^atm ) 

where the superscript “band” denotes that the factor depends on 
the wavelength measured and X is the airmass which is neces¬ 
sary for considering the optical path length of the atmosphere 
during the observation. For the sake of simplicity, we assume 
a homogeneous atmosphere in terms of airmass to obtain the 
atmospheric transmittance at a given airmass by scaling that 
at the airmass of 1.0. In fact, such a simple scaling tends to 
overestimate the atmospheric extinction at X > 1 due to the 
non-linear dependence of the amount of extinction on airmass 
(Manduca & Bell 1979; Tokunaga et al. 2002). Thus, the 
resultant quantities, T atm and PWV, are considered to be lower 
and upper limits, respectively. Tgiter is measured at 77 K in the 
laboratory (Table 2). We calculate using the ATRAN 

model. As the atmospheric transmittances at the H and K a - 
bands are largely independent regardless of PWV (see Figure 
1), the factor x T^nir can t> e calculated, for example. 


for the A' s -band as the following: 


rpK s v rpK s 
^ Tel x 1 ANIR 


K s 

Vobs 


T Ks v (T Ks A 
1 filter * \ atm/ 


K s 


( 2 ) 


In the same way, we calculate the factor for the //-band 
(T^j x Since the factor for the fT-band is much the 

same as that for the A's-band with a dispersion of a few per¬ 
cent, indicating that the factor has small dependence on wave¬ 
length, we then interpolate T^ el x T^ KlR and T .^ x T^ IR 
to estimate the same factor for a NB filter (TV 1875 or TV191), 
x ^anir- Finally, we obtain PWV by iteratively calculat¬ 
ing the left-hand integral in Equation (3) within the bandpass 
[Ai, A 2 ] of the NB filter with changing PWV to make it equal 
to the right-hand value. 


A 2 /T.NB 
JAi "^filter 


/ T PWV, NB 
X V-^atm 


) x dX 





o^NB w t^NB 
^ Tel x 1 ANIR 


( 3 ) 


The evaluated PWV values and the corresponding atmo¬ 
spheric transmittances at the TV 1875 (TV 191) filter are shown 
in the left- and right-hand ordinates of Figure 12a (12b), re¬ 
spectively. Each data point is calculated by using a reduced 
and combined NB image of several dithering frames which has 
a total integration time of ~ 1200 s. The vertical error bars 
include in quadrature (i) the lcr uncertainty of the through¬ 
put of the NB data and (ii) frame-by-frame fluctuations of the 
throughput in the same dithering sequence. We evaluate the 
latter factor by processing individual frames in the same man¬ 
ner as described above for two objects: one with a higher (~ 
0.8 mm) and another with a lower (~ 0.2 mm) PWV evaluated 
from their stacked image. We find a frame-by-frame fluctua¬ 
tion of about 0.14 mm for both objects, and we thus introduce 
a constant factor of 0.14 mm as the frame-by-frame fluctuation 
for all the data in Figure 12. Our result tends to show lower 
PWV values (the median and its lcr dispersion are 0.40 ± 0.30 
for TV1875 and 0.37 ± 0.21 mm for V191) than those reported 
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Fig. 12. Left : Comparison of PWV values measured at the APEX site (5,100 m) with those at the TAO site (5,640 m) evaluated via the atmospheric 
transmittance at the TV 1875 filter. The vertical error bars account for only random errors and are the sum in quadrature of (i) the lcr uncertainties in the 
throughput of the TV 1875 filter which is used to estimate the atmospheric transmittance (the right-hand ordinate) and (ii) the frame-by-frame fluctuation 
of the throughput in the same dithering sequence (0.14 mm in PWV for all the data). The horizontal error bars indicate lcr dispersions of the APEX 
PWV values during TV 1875 observations for each object (typical integration time 1200 s including overheads). The dashed line indicates a one-to-one 
relation and the dotted line shows the best-fitting linear relation with a slope of 0.47. The shaded region shows the possible region expected from the 
exponential distribution of the water vapor with a scale height of 0.3-1.9 km. Right : Same as the left-hand panel, but for the TV 191 data. The best-fitting 
(dotted) line has a slope of 0.60. 


previously (i.e., Giovanelli et al. 2001b; Erasmus & Sarazin 
2002). This may be partly because our observations have been 
carried out mostly in May and October, when the climate is ex¬ 
pected to be drier than the yearly mean. Note that the effect of 
the airmass has been corrected in Figure 12, so that the disper¬ 
sion of each measurement of PWV is likely to be caused by the 
temporal fluctuation in the atmospheric absorption. 

For a quantitative comparison of the result, we use archival 
PWV data measured by a radiometer of Atacama Pathfinder 
Experiment (APEX, Glisten et al. 2006) located near the 
base of Cerro Chajnantor at an altitude of 5,100 m. We extract 
the APEX PWV data during our iV 1875 (1V191) observations, 
which are plotted along the abscissa in Figure 12a (12b). The 
error bars represent 1 a dispersion of the extracted APEX PWV 
during our observations (~ 1200 s per object). We clearly see 
that the PWV evaluated at the TAO site is remarkably lower 
(the median ratio and its la dispersion are 49% ± 38% for 
A71875 and 59% ± 26% for 7V191) than those measured at the 
APEX site. A simple linear fit through the origin to the dis¬ 
tribution yields a slope of 0.47 for 7V1875 and 0.60 for 7V191 
(the dotted line in Figure 12). 

Let us now consider the vertical distribution of the water va¬ 
por to discuss possible causes of the difference in the PWV be¬ 
tween the APEX (PWVapex) and the TAO (PWVtao) sites. 
When the water vapor is assumed to be distributed exponen¬ 
tially, PWV is derived by integrating an exponential profile, 
p = po exp[— {h — ho)/h e ], over altitude h, where p is the 
water vapor density (in kg m 3 ) at a given altitude h above sea 


level (in km), po the density at a reference altitude ho, h e the 
scale height at which the water vapor density decreases by a 
factor of e. Giovanelli et al. (2001b) have measured the verti¬ 
cal distribution of the water vapor above the Chajnantor plateau 
by a combination of radiometric (at both 183 and 225 GHz) and 
radiosonde measurements, and derived the median scale height 
h e ~ 1.13 km by fitting the individual distributions to the ex¬ 
ponential profile with ho of 5.0 km (i.e., Chajnantor Plateau). 
By using a scale height of 1.13 km, the PWV at the TAO site is 
calculated for a given PWV at the APEX site, as follows (see 
also Otarola et al. 2010, 2011). The density po is found by 
integrating the above equation above 5,100 m since the integral 
equals to the PWV at the APEX site. For example, the PWV 
at the APEX site of 1.0 mm leads to po ~ 1.44 g cm -3 . Then, 
the PWV at the TAO site is derived, to be ~ 0.62 mm in the 
case of the example, by integrating the equation above 5,640 m 
with po substituted. Using the measurements of the exponen¬ 
tial scale height at night time (the peak-to-peak values in h e 
of ~ 0.3-1.9 km) by Giovanelli et al. (2001b), we derive the 
corresponding PWV values at the TAO site as a function of the 
PWV at the APEX site, which are shown in Figure 12 with a 
shaded region. We find that almost all of the data points are in 
excellent agreement, within the uncertainties, with the expecta¬ 
tion of such exponential distributions of the water vapor having 
the measured scale heights. Note that there are a few outliers 
above the region, which might be caused by liquid phase of 
water (fog or clouds) in the atmosphere which is not detected 
by the APEX radiometer operating at 183 GHz, as suggested 
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by previous studies (Matsushita & Matsuo 2003; Tamura et al. 
2011). We should also note the possible influence of the pres¬ 
ence of temperature inversion layers on the water vapor distri¬ 
bution. Giovanelli et al. (2001b) find that their radiosonde data 
often show temperature inversions which make the distribution 
far from the exponential shape, so that much of the water vapor 
would be trapped below the inversion layers. To explore the 
influence requires a large sample of radiosonde data (or any 
equivalent data on the vertical distributions of the water vapor) 
and it is far beyond the scope of this paper, but we consider 
that our conclusion of the advantage of the site remains unre¬ 
vised, because temperature inversions often take place below 
the altitude of 5,500 m at night time (Giovanelli et al. 2001b), 
indicating systematically lower PWV content at the TAO site 
than at the APEX site, which is qualitatively the same trend as 
derived from Figure 12. 

While there is relatively large uncertainties in our analysis, it 
is encouraging to see that the low PWV content at the TAO site 
is confirmed by independent methods (i.e., NB imaging and 
radiosonde), which suggests that the TAO site at the summit 
of Cerro Chajnantor is suitable for infrared astronomy, and in 
particular that miniTAO/ANIR has an excellent capability for 
Paa observations at around A = 1.9 //m. 

5. Summary 

We have developed a near-infrared camera, ANIR (Atacama 
NIR camera), for the 1.0-m miniTAO telescope installed at the 
summit of Cerro Chajnantor (an altitude of 5,640 m) in north¬ 
ern Chile. The camera covers a field of view of 5(1 x 5(1 with 
a spatial resolution of 0"298 pixel 1 in the wavelength range 
of 0.95 to 2.4 /itn. 

The unique feature of the camera coupled with advantages 
of the site is a capability of narrow-band imaging observations 
of a strong hydrogen emission line, Paschen-a (Paa), at A = 
1.8751 pm from the ground, at which wavelength it has been 
quite difficult to conduct ground-based observations so far due 
to deep atmospheric absorption mostly from water vapor. We 
have been successfully obtaining Paa images of Galactic ob¬ 
jects and nearby galaxies since the first light observation with 
ANIR in 2009. 

The throughputs at the narrow-band filters (iVl875 and 
TV 191) show larger dispersion (~ 10%) than those of the broad¬ 
band filters (a few percent), indicating that they are affected 
by temporal fluctuations in Precipitable Water Vapor (PWV) 
above the site. Combining the atmospheric transmission model 
with the throughputs, we evaluated the atmospheric transmit¬ 
tance at the narrow-band filters and the PWV content at the 
site. We find that the median and the dispersion of the PWV 
are 0.40 ± 0.30 for iV1875 and 0.37 ± 0.21 mm for AT191 
data. Comparing those data with the radiometer data taken 
by APEX, we find that the site has remarkably lower (49% ± 
38% for A/1875 and 59% ± 26% for N 191) PWV values than 
the APEX site (5,100 m). The differences in PWV between 
those sites are found to be in excellent agreement with those 
expected from the exponential distribution of the water vapor 
with scale height within 0.3-1.9 km which has been measured 
using the radiosonde at night time (Giovanelli et al. 2001b). 
Although temperature inversions (which mostly take place be¬ 


low the summit at night-time) make the water vapor distribu¬ 
tion far from exponential, it leads to lower PWV content above 
the inversion layers, which is qualitatively consistent with our 
findings described above. Taken all together, we conclude that 
miniTAO/ANIR and the site, the summit of Cerro Chajnantor, 
provides us with an excellent capability for ground-based Paa 
observations. 
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